Spatial Landau-Zener-Stiickelberg interference in spinor Bose-Einstein condensates 
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We investigate the Stiickelberg oscillations of a spin-1 Bose-Einstein condensate subject to a 
spatially inhomogeneous transverse magnetic field and a periodic longitudinal field. We show that 
the time-domain Stiickelberg oscillations result in modulations in the density profiles of all spin 
components due to the spatial inhomogeneity of the transverse field. This phenomenon represents 
the Landau-Zener-Stiickelberg interference in the space-domain. Since the magnetic dipole-dipole 
interaction between spin-1 atoms induces an inhomogeneous effective magnetic field, interference 
fringes also appear if a dipolar spinor condensate is driven periodically. We also point out some 
potential applications of this spatial Landau-Zener-Stiikelberg interference. 

PACS numbers: 03.75.Lm, 03.75.Mn, 67.85.Fg 
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I. INTRODUCTION 

Quantum two-level systems often exhibit an avoided 
energy-level crossing which can be traversed using an 
external control parameter. If one sweeps the exter- 
nal parameter through the avoided crossing, a coherent 
Landau-Zener transition occurs [H, Hj]. When traversing 
the avoided crossing twice, by sweeping the external pa- 
rameter back, the dynamical phase accumulated between 
transitions may give rise to constructive or destructive 
interference in the time-domain, known as Stiickelberg 
oscillations [!, 0] ■ When such a two- level system is sub- 
jected to a periodic driving in time, the physical ob- 
servables of the system exhibit a periodic dependence 
on some external parameters, which is referred to as 
Landau-Zener-Stiickelberg interferometry Recently, 
it was shown that Landau-Zener-Stiickelberg interfer- 
ometry is of particular importance to superconducting 
qubits la-tlCl. nitrogen vacancy center [Tl|, and quan- 
tum dots ll3T | , as it provides an alternative means 
to manipulate and characterize the structure of two-level 
systems. 

In atomic physics, Stiickelberg oscillations were ob- 
served in the dipole-dipole interaction between Ryd- 
berg atoms with an externally applied radio-frequency 
field [3]. In particular, for ultracold atomic gases, 
Stiickelberg oscillations were demonstrated using, the 
internal-state structure of Feshbach molecules [HI, [l6[ 
and using Bose-Einstein condensates in accelerated op- 
tical lattices 17| . There are numerous theoretical works 



studying Landau-Zener tunneling subject to a temporal 
periodic driving for condensates trapped in double-well 
potentials |18l - [2fj |. Vasile et al. [21j also proposed an 
interferometer of spinor condensates using Stiickelberg 
oscillations. 

In this work, we study the dynamics of a spin-1 conden- 
sate subject to a spatially inhomogeneous transverse mag- 
netic field. When the condensate is driven by a tempo- 
rally periodic longitudinal magnetic field, the Stiickelberg 
phases accumulated at different spatial positions are dif- 



ferent, which results in modulations in the density pro- 
files of all spin components. Therefore, by imposing a 
spatially non-uniform transverse field, we convert the 
time-domain interference into a space-domain one, which 
we refer to as spatial Landau-Zener-Stiickelberg interfer- 
ence (SLZSI). In spinor condensates, the spatially inho- 
mogeneous transverse field can also be provided by the 
magnetic dipole-dipole interaction between atoms. We 
show that the SLZSI occurs even in the absence of any 
external transverse field. This phenomenon can be used 
to detect dipolar effects in a spinor condensate. 



For a condensate confined in a Ioffe-Pritchard trap, 
Leanhardt et al. [22j experimentally demonstrated that 
a vortex can be imprinted by adiabatically inverting the 
axial magnetic field. Interestingly, they also swept the 
axial field back to its original direction. Due to the adia- 
baticity of the whole process, they found that the conden- 
sate recovered its original state. The dynamics of a spin-1 
gas subject to a temporally oscillating field was studied 
in various papers |23l - l25| . However, to the best of our 
knowledge, we are not aware of any work on Stiickelberg 
oscillations subject to a spatially inhomogeneous mag- 
netic field. 



This paper is organized as follows. In Sec. [TT1 we 
consider a continuum of two-level systems subject to a 
non-uniform transverse field and a periodic driving along 
the longitudinal direction. We show that the occupation 
probabilities are periodically modulated in the space- 
domain. Section Hill studies the SLZSI of a spin-1 conden- 
sate in a Ioffe-Pritchard trap. We show that the positions 
of the destructive interference points agree with those 
presented in Sec. |TTJ In Sec. |TV] we study the SLZSI in- 
duced by the magnetic dipole-dipole interaction. Finally, 
we present our conclusions in Sec. [V] 
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FIG. 1: (Color online) Schematic plot of the adiabatic en- 
ergy levels of a spin-i system subject to a transverse field 
and a temporally periodic driving field along the longitudinal 
direction. 



II. SPATIAL QUANTUM INTERFEROMETRY 
IN THE SINGLE-PARTICLE PICTURE 

Let us first briefly summarize the Landau- Zener tun- 
neling of a single two-level system under temporally peri- 
odic driving [5] . This will allow us to introduce quantities 
that will be used afterwards. Specifically, we consider the 
model Hamiltonian 



H^(t)/tl=-^a x - £ M az 



(1) 



where HA is a constant representing the level splitting, 



e(t) = Acosflt 



(2) 



is the periodic driving with amplitude A and frequency 
f2, and a x ^ z are the Pauli operators. For a spin-i atom, 
A and e(t) can be induced by a constant transverse 
magnetic field and an ac longitudinal field, respectively. 
Throughout this work, we will assume that the longitu- 
dinal periodic driving contains no dc component. 

For the model Hamiltonian Eq. ([IJ, we are interested 
in the time-dependence of the occupation probabilities in 
the upper and lower energy levels at time t. Although 
this problem can be solved by numerous approaches, 
here we quote the results from the adiabatic-impulse 
model (see, e.g., Ref. Q). In Fig. [1] we schematically 
plot the adiabatic energy levels of the system; namely, 
the instantaneous eigenvalues of the Hamiltonian ([1}. 
The avoided-level crossings (with energy splitting HA) 
at times t = (n + i) n/fl are induced by the transverse 
field, where n = 0, 1, 2, • • • . The evolution of the system 
can be divided into two stages: 

• Away from the avoided crossings, the system 
evolves adiabatically by following the evolution ma- 



trix 



where 



U = 



e -«/a o 
e<l 2 



C = - jf ^A* + e{t)Ht 



(3) 



and C/2 is the dynamical phase. 

• Landau-Zener transitions occur at the avoided 
crossings, which can be represented by the tran- 
sition matrix 



-Vp 

VT=pe l<p 



(4) 



where the transition probability is p = exp(— 2tt<5) , 
with (5 = A 2 /(4A£l), and the phase jump is tp = 
- f + 5 (In S - 1) + argr(l - i5) , with T(- ) being the 
gamma function. In particular, in the fast-passage 
limit (5 <C 1), the dynamical phase and phase jump 
become, respectively, £ « A/£l and ip « — 7r/4. 

After the two-level system is driven for £ half-periods, 
one may take the measurement of the population ei- 
ther at time t\ or t 2 as shown in Fig. [1] Correspond- 
ingly, the total evolution matrix becomes Ui(UNU) or 
U 2 U~ 1 {UNU) l+1 , respectively. One immediately sees 
that, for the purpose of calculating the transition prob- 
ability, a very relevant quantity is the transition ma- 
trix m, m 



(UNU) 1 = Un "" 21 

1 "21 U* n 



where 



W 1 -psm(</j + C) 



sin(^0) 



sin f 



"21 

cos 



Vp 



sm{lO) 



sin ( 



9 = V^ ~pcos(ip + (). 



(5) 

(6) 

(7) 
(8) 



Assuming that only the upper level is initially populated, 
the occupation probabilities at the end of the field sweep- 
ing become 



P+ = cos 2 (£e) + (1 -p)sin% + C)— ., 

sin 

p = 

sin 2 9 



(9) 
(10) 



Now we generalize the above single two-level system 
to a continuum of isolated two-level systems which are 
distributed over the range x € [0, L\. We further assume 
that the energy splitting depends linearly on the position 
x, i.e., 



A(x) = b'x + b , 



(11) 
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FIG. 2: (Color online) (a) Schematic plot of the right-hand- 
side (solid line) and the left-hand-side (horizontal dashed 
lines, corresponding to different fc's for a given q) of Eq. (|12|l . 
Destructive inferences occur at the Xi's. (b) Spatial distribu- 
tion of the transition probability P_ [Eq. (|10[) ] for the follow- 
ing dimensionless parameters: b' — 0.01A/L, bo = 0, fl — A, 
and q = 100. 



where b' and bo are two constants. Without loss of gen- 
erality, we assume that b' > and b > 0. Due to the 
position dependence of the energy splitting, the transi- 
tion probabilities now must depend on the position of the 
two-level system. To proceed further, let us focus on the 
occupation probability P_. If P_ is measured after the 
periodic driving is applied for q periods (i.e., i = 2q), 
destructive interference for P_ occurs at 2q6{x) = kn 7 
with k being an integer. Using Eq. (|8]l. the condition for 
destructive interference becomes 



kn 



cos — = y/l — p(x) cos[ip(x) + ((x)] 
2q 



(12) 



In Fig. [2ja), we schematically plot the left- and right- 
hand-side of Eq. (|12p . The x coordinates of the inter- 
sections of the left- and right-hand-side then represent 
the positions of different destructive interferences. As 
a consequence, non-trivial spatial structure forms in the 
transition probability, representing the SLZSI. 

To gain more insight into the spatial structure of the 
transition probability, we consider the fast-passage limit 
by assuming that 5(x) <C 1 for x G [0, L]. The destructive 
interference condition can be approximated as 



kn 

cos — ~ CA{x), 
Iq 



(13) 



where C = yj-rr/ (2AQ) cos (— f + n) 1S a constant. In 
the fast-passage regime, we may assume that A(x) is 
sufficiently small, which allows us to focus on the val- 
ues of k in the vicinity of q. To this end, we rewrite k as 



--q + k', with k' = 0,±1,±2,- • • 
<C q, equation (ITBl reduces to 

k'ir 



By further assuming 



— ~C(b'x + b ). 
2q 



(14) 



Apparently, corresponding to different k' , the transition 
probability are spatially modulated with equal interval 



(15) 



2qV\C[ 

In addition, the necessary condition for the range L to 
accommodate a destructive interference is L > l x , which 
implies that the temporal periodic driving must be ap- 
plied for over ir/(2b'L\C\) periods. 

In Fig. [2fb), we illustrate an example of the spatial 
distribution of P_, which is plotted using Eq. (fl"U|) for 
the set of parameters b' = 0.01^4/P, 6o = 0, ft = A, and 
q = 100. As can be seen there, P~(x) exhibits nice spatial 
periodicity. The spatial period read out from Fig. [^b) is 
in very good agreement with that predicted by Eq. f| 15[) . 
The nearly sinusoidal dependence of P_ on x can be in- 
tuitively understood as follows. In the limit 8(x) -C 1, p 
(w 1) and 8 have a very weak position-dependence such 
that sin# is essentially a constant. However, when I is 
sufficiently large, even a weak x-dependence in 6 is sig- 
nificantly amplified in the function sin 2 (£9), which is the 
only term in P_ contributing to the position dependence. 



III. SLZSI OF A SPIN-1 CONDENSATE IN A 
MAGNETIC TRAP 

In this section, we study the SLZSI of an optically 
trapped spin-1 condensate subject to a transverse mag- 
netic field of the form of a Ioffe-Pritchard trap and a 
temporally periodic driving along the longitudinal direc- 
tion. To this end, we first consider the spin dynamics of a 
single atom in the F = 1 hyperfine state subject to a con- 
stant transverse magnetic field B x and a time-dependent 
longitudinal field B z (t) = Bq cos tit. The Hamiltonian of 
the system reads 



H^(t)/fi=-AF x -e(t)F z , 



(16) 



where A = —gp[iBB x jti and e(t) = —gp^iBBoCOsQt/h, 
with i?f(= — 1/2) being the Lande g-factor of the atom, 
fiB the Bohr magneton, and F = (F x , F y ,F z ) the angular 
momentum operator. Note that a spin-1 atom contains 
three Zeeman sublevels, corresponding to the magnetic 
quantum number a — 1,0, and —1. Utilizing the Ma- 
jorana representation 0, |27I |. the spin dynamics of the 
Hamiltonian (| L6[) can be easily derived from that of a 
spin-i particle. Specifically, the evolution matrix after I 
half-periods becomes 



y/2- 



'ii 

U11U21 



— V2uilU2i 

\un\ 2 - \u 2 i\ 2 
V2unU 2 i 



-V2u 

Ki) 



ll u 21 



(17) 
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where itn and U21 are given, respectively, in Eqs ([6]) and 
([7]). For an initial state with only the a — 1 spin com- 
ponent populated, the occupation probabilities then be- 
come Pi = Ki| 4 , P = 2\un\ 2 \u 2 i\ 2 , and P_i = |u 2 i| 4 - 
Obviously, there exists a relation 

P = 2^P 1 P- 1 (18) 

among the occupation probabilities of different spin com- 
ponents. 

For an ultracold atomic gas, atoms interact with each 
other and are free to move. Therefore, to strictly realize 
the SLZSI obtained from isolated two-level systems, the 
effects of the kinetic and interaction energies have to be 
eliminated. To achieve this, one may load the atoms into 
an optical lattice potential. If the depth of the optical po- 
tential is so high that the system is in the Mott insulator 
phase, the kinetic energy can be safely ignored. Further- 
more, in the case where each lattice site is only occupied 
by a single atom, the atom-atom interaction can also be 
eliminated. However, as we shall show below, for typi- 
cal experimental parameters, even for a singly trapped 
spinor condensate, SLZSI can be well described by the 
model introduced in Sec. ITU 

Now we proceed to study the SLZSI of a spin-1 conden- 
sate. Atoms interact with each other via the short-range 
potential [H, [H 

V (r - r') = (c + c 2 F • F')5(r - r'), (19) 

where cq = 47rfi, 2 (ao + 2a2)/(3m) and C2 = 47rS 2 (a2 — 
ao)/(3m), with a/ (/ = 0,2) being the s-wave scatter- 
ing length in the combined symmetric channel of total 
spin /. For the sodium atoms considered in this sec- 
tion, ao = 50ob and a 2 = 55as, with as being the Bohr 
radius. Note that the parameters Co and C2 represent, 
respectively, the strength of the spin-independent and 
spin-exchange collisional interactions. The spin-exchange 
interaction of the sodium atoms is antifcrromagnetic. 
There also exist atoms whose spin-exchange interaction 
is ferromagnetic (e.g., rubidium). However, for the mag- 
netic field considered in this section, the Zeeman energy 
is much larger than the spin-exchange interaction energy 
such that the spin-exchange interaction is unimportant 
to the spin dynamics. 

We assume that the transverse magnetic field takes the 
form of a Ioffe-Pritchard trap, such that the total external 
field becomes 

B ext {r,t) = B'(xk-yy) + B z (t)z, (20) 

where B 1 is the gradient of the transverse magnetic field. 
Furthermore, to stably confine the condensate, a spin- 
independent optical trap 

^opt(r) = ^muj'i{x 2 + y 2 + A 2 z 2 ) 

is also applied, where u>± is the radial trap frequency and 
A is the trap aspect ratio. For the numerical simulations 



presented below, we shall choose B' = 15G/cm, Bq = 
1G, ui x = (2tt)100Hz, and A = 6. 

Within the framework of mean- field theory, a spin-1 
condensate of N atoms is described by the wave functions 
^ Q (r) (a — 0,±1), which satisfy the dynamic equations 

ih-Qf = [T+ U opt + c n(r)]ip a + g F fi B B cS ■ ¥ a pxj)p, 

(21) 

where T = — ft 2 V 2 /(2m) is the kinetic energy term, 
n(r) = ip^ipa is the number density of the condensate 
normalized to the total number of atoms N, and 

B cff (r, t) = B«t(r, t) + B cxc (r, t) (22) 

is the effective magnetic field which includes the external 
magnetic field and the contribution originating from the 
spin-exchange interaction 

B« c (r,t) = -^-S(r,t) (23) 

9FPB 

with S(r,i) = Y^ap^V^apipp being the spin density. 
Here we have neglected the magnetic dipole-dipole inter- 
action since it is much smaller than the Zeeman energy in 
a Ioffe-Pritchard trap. However, the effect of the dipolar 
interaction will be addressed in the next section. 

We demonstrate the SLZSI in a magnetic trap by nu- 
merically evolving Eqs. (l2lT) for a condensate of N = 
2 x 10 6 sodium atoms. The initial wave functions are 
taken as the ground state of the spinor condensate under 
the external field B ext (r, 0). In the results presented be- 
low, we will focus on the behavior of the column density 
of the different spin components 

n a (x,y) = J dz\^ a {r)\ 2 , (24) 

which also corresponds to the absorption image of the 
atomic gas. 

Figure [3] shows the column densities of all spin com- 
ponents for Q = 107ru;j^ and after the driving field being 
applied for q = 0, 2, 4, and 6 periods. Due to the ax- 
ial symmetry of the column densities, they are plotted 
as functions of p — y/ ' x 2 + y 2 . Initially, only the a = 1 
component is populated for the given parameters. As 
we start to drive the condensate with an ac longitudinal 
field, other spin states also become occupied. In partic- 
ular, ripples start to develop in the density profiles of all 
spin components. If the condensate is driven for a longer 
time, more ripples will appear. 

The positions of the destructive interference in the 
a = — 1 component is also determined by Eq. ([T2|) . Here, 
instead, we will focus on fi\ since the a = 1 component 
contains the majority of the atoms when the system is 
only driven for a short period of time. For the given 
parameters, we have 1 — p(p) <C 1. Therefore, the po- 
sitions of the destructive interference in n\ are mainly 
determined by the first term on the right-hand-side of 
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FIG. 3: (Color online) Column densities at various times for 
Q, — 1Qhui±. The solid, dashed, and dash-dotted lines de- 
note, respectively, the a = 1, 0, and —1 spin components. 
The vertical arrows denote the positions of the destructive 
interferences predicted by Eq. (|25l) . 



Eq. (|9|). Following the same analysis as that presented 
in Sec. HI1 destructive interference in n\ occurs at the 
positions determined by the equation 

2k' - 1 



= y/l-p(p)cos[<p(p) + C(p)] (25) 

with k' being integer, where q is the number of periods 
that the driving field has been applied for. For the pa- 
rameters given in Fig. [3l the location of the destructive 
interference points predicted by Eq. (|2"5j) are shown as 
black vertical arrows in the figure. As can be seen, they 
agree well with those obtained through the full numerical 
simulations. 

To gain more insight into the SLZSI of a spin-1 con- 
densate, let us verify the relation Eq. (TT51) . which holds 
rigorously in the single-particle picture. For the col- 
umn densities of the condensates, Eq. (|18|) becomes 
no = 2^/n\n~[. In our numerical simulation, we find 
that tt-o and 2y / nTn3j" agree with each other for small t. 
However, as shown in Fig.Uta), a prominent discrepancy 
appears for uj±t = 1. Since Eq. (fl"8)l is derived from the 
Majorana representation, its violation indicates a devia- 
tion from the single-particle picture. To identify which 
source causes the discrepancy, one may evolve Eq. (|21[) 



c 
is 





(a) 


1 /ill 
1 1 If 1 

1 K /, I 

1 VI 1 

/ Xk 


1 V 



10 20 

p(pra) 






30 



FIG. 4: (Color online) Verification of the relation Eq. (fl8|) 
with (a) and without (b) the contribution from the kinetic 
energy for ui±t — 1. The solid and dashed lines correspond 
to 2-y/nin_i and no, respectively. Other parameters are the 
same as those used in Fig. [3] 



with the kinetic energy term removed. This treatment is 
equivalent to taking the Thomas-Fermi approximation. 
The result presented in Fig. @Jb) shows that 2- s /nin_i is 
visually identical to uq after the kinetic energy term is 
dropped, which suggests that the discrepancy in Fig.Hta) 
is mainly caused by the center of mass motion of the 
atoms. 

Vengalattore et al. [30] have demonstrated that spinor 
condensates can be regarded as high-resolution magne- 
tometers. Potentially, due to its sensitivity to the inho- 
mogeneity of the magnetic field, the SLZSI in a spinor 
condensate can also be used to measure the gradient of 
the magnetic field. 



IV. SLZSI IN A DIPOLAR SPIN-1 
CONDENSATE 

Now, we turn to study the SLZSI in a dipolar spinor 
condensate. In addition to the short-range interaction 
V , there also exists a long-range magnetic dipole-dipole 
interaction between two spin-1 atoms. The interaction 
potential between two magnetic dipoles takes the form 



V d (r-r') = c d 



F-F'-3(F-e)(F'-e) 



r !\3 



(26) 



where the strength of the dipolar interaction is charac- 
terized by Cd = pop%gp/(4:Tr), with po being the vacuum 
magnetic permeability and e = (r — r')/|r — r'| is a unit 
vector. Within the framework of mean-field theory, the 
dipolar interaction generates an effective magnetic field 
of the form [3ll 



BdipM) 



dr 



,S(r ',t)-3[S(r',t) 

It- - t-'|3 



(27) 



9FPB 

Consequently, the total effective magnetic field becomes 
B off (r, t) = B ext (r, t) + B oxc (r, t) + B dip (r, t). (28) 
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FIG. 5: (Color online), (a) Magnitude of the transverse field 



B 



dtp O 11 un its of fiG) induced by the dipolar interaction for 

a spin polarized condensate, (b) p-dependence of B^ for 
various z = zq planes. 



We note that the magnitude of the transverse component 
of Bdip can be formally expressed as 



D 



(-L) = 
dip 



B2 ip -(B dip .z)2 



(29) 



As we will show, B^j is non-uniform. Therefore, SLZSI 
may be induced in spinor condensates even in the absence 
of an external transverse magnetic field. 

To proceed further, we consider a concrete example of 
a spin-1 condensate containing N = 10 7 rubidium atoms. 
The Lande g- factor of the atoms is gF = —1/2. Moreover, 
the s-wave scattering lengths between rubidium atoms 
are a>Q = 101 .8a b and a 2 = 100.4as. The optical trap 
has the same parameters as those adopted in Sec. Mil 
Finally, to emphasize the effect of the dipolar interac- 
tion, we assume that the external field only contains a 
longitudinal component, i.e., 

Bext(t) = Bq cosfiiz. 

The value of Bq is typically around several hundreds pG, 
under which the spins of the atoms are fully polarized 
along the z-axis. 

Before we turn to study the dynamics of the conden- 
sate under an external driving field, it is instructive to 
examine the structure of the magnetic field induced by 
the dipolar interaction. For simplicity, we assume that 
only the a = 1 spin component is populated, which 
is essentially the ground state under the external field 
B ex t(^ = 0). Figure a) shows the magnitude of the 
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FIG. 6: (Color online) Column densities at various time t's 
for Bo = 100 /iG and Q = 0.2nuj±. The solid, dashed, and 
dash-dotted lines denote, respectively, the a = 1, 0, and — 1 
spin components. The dotted lines represent the total column 
densities. 



transverse component of Bdi p for a spin-polarized con- 
densate. Due to the cylindrical symmetry of the system, 
fl^ip reduces to a function of p and z. Clearly, B^j is 
non-uniform and takes a butterfly shape in the pz-plane. 
For this specific example, the maximum value of B^j is 
around 10.4 fiG. In particular, as can be deduced from 
Eq. (|2T|) . the effective transverse field vanishes in the xy- 
plane. 

To reveal more details about the transverse field, we 
plot B^rj(p,z — z Q ) for various z 's in Fig. [DJb). On 

a given z — zq plane, B^(p, zo) is roughly a linear 
function when p is small. However, the gradient of the 
transverse field sensitively depends on the value of zq . 
Moreover, B^ becomes a time-dependent function once 
the external driving field B cxt is applied. Therefore, one 
should not use Eq. (f25|) to predict the positions of the 
destructive-interference points. 

The dynamics of the condensate can be simulated by 
numerically evolving Eqs. (1211) with the effective mag- 
netic field in Eq. (|28[) . The initial wave function is taken 
as the ground state under the external field B cxt (i = 0). 
In Fig.O we present the column densities of all spin com- 
ponents after the driving field is applied for various peri- 
ods. The parameters of the driving field are Bq — 100 pG 
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FIG. 7: (Color online) Column densities at time uj±t = 200 
for B = 500 fiG with SI = 0.2tkj ± (a) and SI = 0.04ttwx (b). 
The solid, dashed, and dash-dotted lines denote, respectively, 
the a = 1, 0, and —1 spin components. The dotted lines 
represent the total column densities. 

and SI — 0.2ttll!±. As can be seen, even though the to- 
tal density remains unchanged, ripples start to develop 
on the density profiles of all spin components after the 
driving field is applied for a few periods, which indicates 
SLZSI occurs in a spin-1 condensate subject to a peri- 
odic driving. More ripples will appear if one evolves the 
system for a longer time. 

The structure of the density ripples also depends sen- 
sitively on the parameters of the driving field. As it can 
be seen from Fig. EJa), the ripples in hi are significantly 
suppressed if we increase the amplitude of the driving 
field to B = 500 fiG. This can be intuitively understood 
as follows. A destructive interference on hi occurs when 
a considerable number of atoms in the a = 1 compo- 
nent are transferred to other spin components. If one 
increases Bq while keeping SI unchanged, the transition 
probability decreases. Therefore, in order to gain visi- 
ble density ripples, one has to lower the frequency of the 
driving field. Indeed, as shown in Fig. [T^b) , density rip- 
ples appear again by lowering the driving frequency to 
SI = 0.04tto; ± . 

From Eq. (j2"3")l . it is tempting to think that, even in 
the absence of the magnetic dipole-dipole interaction, the 
spin-exchange interaction can also induce density ripples. 
This turns out to be untrue for the initial states consid- 
ered in this section. To show this, we consider an initial 
state with spins being polarized to an arbitrary direc- 
tion. For such a state, the total spin of the condensate is 
N and the number of atoms in all spin components are 
uniquely determined [32l. l33j 

N ( M\ 2 , N ( M 2 \ 

where M — Ni — N^i is the ^-component of the to- 
tal spin. Now we assume that an magnetic field is ap- 
plied along the z-axis. Since the total spin (N) and its 



^-component (M) are conserved in the absence of the 
dipolar interaction |28j], the contact spin-exchange inter- 
action will not cause any spin-mixing. Therefore, for an 
initially polarized spin-1 condensate subject to a peri- 
odic driving, the appearance of density ripples is an un- 
ambiguous manifestation of the dipolar interaction in a 
spinor condensate. 

In a previous work [iU , we proposed to detect the effect 
of dipolar interaction in spinor condensates by adiabati- 
cally inverting the longitudinal field, where the strength 
of the longitudinal field is only around several tens of 
micro-Gauss. Compared to that scheme, the obvious ad- 
vantage of using the dipolar-interaction-induced SLZSI is 
that the strength of the longitudinal field can be much 
higher. 



V. CONCLUSION 

We have proposed a method to convert the time- 
domain Stiickelberg oscillations to an interference pat- 
tern in the space-domain by imposing a spatially non- 
uniform transverse magnetic field. For a continuum of 
two-level systems, we showed that the occupation prob- 
abilities are periodically modulated in space. In addi- 
tion, we also obtained a relation between the spatial pe- 
riod and the system parameters in the fast-passage limit. 
We then demonstrated the SLZSI for a spin-1 conden- 
sate subject to a transverse field of the form of a Ioffe- 
Pritchard trap. We found that the kinetic and interaction 
energies only slightly modify the interference patterns ob- 
tained from the single-atom model if the system is not 
driven for a long time. Finally, we showed that the SLZSI 
can also be induced by the magnetic dipole-dipole inter- 
action, even in the absence of an external transverse field. 
Potential applications of the SLZSI include the measure- 
ment of the spatial inhomogeneity of the magnetic field 
and the detection of weak magnetic dipole-dipole inter- 
actions in spinor condensates. Finally, we want to point 
out that SLZSI is a single-particle property, therefore, 
a fermionic gas should also exhibit the phenomena de- 
scribed in this work. 



Acknowledgments 

We thank S. Ashhab for valuable comments on the 
manuscript. FN acknowledges partial support from 
the Laboratory of Physical Sciences, National Security 
Agency, Army Research Office, National Science Founda- 
tion under Grant No. 0726909, DARPA, AFOSR, JSPS- 
RFBR under Contract No. 09-02-92114, Grant-in-Aid 
for Scientific Research (S), MEXT Kakenhi on Quantum 
Cybernetics, and the Funding Program for Innovative 
R&D on S&T (FIRST). CPS acknowledges the supports 
from the NSFC under Grant Nos 10935010. SY acknowl- 
edges the supports from the NSFC (Grant Nos. 11025421 



8 



and 10974209) and the "Bairen" program of the Chinese Academy of Sciences. 



L. Landau, Phys. Z. Sowjetunion 2, 46 (1932). [18 

C. Zener, Proc. R. Soc. London, Ser. A 137, 696 (1932). 
E. C. G. Stiickelberg, Helv. Phys. Acta 5, 369 (1932). [19 
E. Majorana, Nuovo Cimento 9, 43 (1932). 
S. N. Shevchenko, S. Ashhab, and F. Nori, Phys. Rep. [20 
492, 1 (2010). 

W. D. Oliver, Y. Yu, J. C. Lee, K. K. Berggren, L. S. [21 
Levitov, T. P. Orlando, Science 310, 1653 (2005). 
M. Sillanpaa, T. Lehtinen, A. Paila, Y. Makhlin, and P. [22 
Hakonen, Phys. Rev. Lett. 96, 187002 (2006). 
S. Ashhab, J. R. Johansson, and F. Nori, Phys. Rev. A 
74, 052330 (2006). [23 

D. M. Berns, M. S. Rudner, S. O. Valenzuela, K. K. 
Berggren, W. D. Oliver, L. S. Levitov, and T. P. Orlando, [24 
Nature 455, 51 (2008). [25 
G. Sun, X. Wen, B. Mao, J. Chen, Y. Yu, P. Wu, S. Han, 

Nat. Commun. 1:51 DOI:10.1038/ncommsl050 (2010). [26 

G. D. Fuchs, V. V. Dobrovitski, D. M. Toyli, F. J. Here- 
mans, and D. D. Awschalom, Science 326, 1520 (2009). [27 

H. Ribeiro and G. Burkard, Phys. Rev. Lett. 102, 216802 [28 

(2009) . [29 
J. R. Petta, H. Lu, and A. C. Gossard, Science 327, 669 

(2010) . [30 
C. S. E. van Ditzhuijzen, A. Tauschinsky, and H. B. van 
Linden van den Heuvell, Phys. Rev. A 80, 063407 (2009). 

M. Mark, T. Kraemer, P. Waldburger, J. Herbig, C. Chin, [31 
H.-C. NSgerl, and R. Grimm, Phys. Rev. Lett. 99, 113201 
(2007). [32 
M. Mark, F. Ferlaino, S. Knoop, J. G. Danzl, T. Kraemer, 
C. Chin, H.-C. N&gerl, and R. Grimm, Phys. Rev. A 76, [33 
042514 (2007). 

A. Zenesini, D. Ciampini, O. Morsch, and E. Arimondo, 
larXiv:1010.2413 (2010). 



G. -F. Wang, L.-B. Fu, and J. Liu, Phys. Rev. A 73, 
013619 (2006). 

X. Luo, Q. Xie, and B. Wu, Phys. Rev. A 77, 053601 
(2008). 

Q. Zhang, P. Hanggi, and J. Gong, Phys. Rev. A 77, 
053607 (2008). 

R. Vasile, H. Makela, and K.-A. Suominen, 
larXiv:0812.0499l (2008). 

A. E. Leanhardt, A. Gorlitz, A. P. Chikkatur, D. Kielpin- 
ski, Y. Shin, D. E. Pritchard, and W. Ketterle, Phys. 
Rev. Lett. 89, 190403 (2002). 

H. Pu, S. Raghavan, and N.P. Bigelow, Phys. Rev. A 61, 
023602 (2000). 

B. Sun and L. You, Phys. Rev. Lett. 99, 150402 (2007). 
K. Gawryluk, K. Bongs, and M. Brewczyk, 
larXiv:1010.0165l (2010). 

S. Ashhab, J. R. Johansson, A. M. Zagoskin, F. Nori, 
Phys. Rev. A 75, 063414 (2007). 

F. Bloch and I. I. Rabi, Rev. Mod. Phys. 17, 237 (1945). 

T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 

(1998). 

M. Vengalattore, J. M. Higbie, S. R. Leslie, J. Guzman, 
L. E. Sadler, and D. M. Stamper-Kurn, Phys. Rev. Lett. 
98, 200801 (2007). 

J.-N. Zhang, L. He, H. Pu, C.-P. Sun, and S. Yi, Phys. 
Rev. A 79, 033615 (2009). 

H. Pu, C. K. Law, S. Raghavan, J. H. Eberly, and N. P. 
Bigelow, Phys. Rev. A 60, 1463 (1999). 
S. Yi, 6. E. Miistecaphoglu, C.-P. Sun, and L. You, Phys. 
Rev. A 66, 011601 (2002). 



